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Abstract, 


A control  system  was  designed  to  improve  the  heading  response  of 
the  Air  Cushion  Landing  System  aircraft  to  differential  thrust  input. 

The  pilot  can  then  control  the  aircraft  in  a more  effective  manner 
while  taxiing  in  gusting  crosswinds.  The  aircraft  equations  of  motion 
and  low  speed  stability  derivatives  were  obtained  from  work  done  at 
the  Air  Force  Flight  Dynamics  Laboratory  at  Wrlght-Patterson  Air  Force 
Base,  Ohio.  Five  different  values  of  engine  response  time-delay  were 
selected  for  the  design  study.  These  time-delay  values  were  0.2,  0.4, 
0.6,  0.3,  and  1.0  seconds,  The  principle  used  for  the  design  was  to 
compenoate  the  system  ignoring  the  time-delay  and  then  to  compensate 
for  the  time-delay.  A two  stage  lead  network  was  used  to  reduce  the 
effect  of  the  engine  time-delay.  The  system  was  simulated  by  using 
a d'.git.al  computer  program  called  MIMIC,  which  simulates  the  functions 
of  an  analog  computer.  The  average  mean  square  error  waB  computed  for 
the  system  with  and  without  the  control  installed,  for  each  of  the 
five  time-delay  values.  In  all  caseB,  the  average  mean  square  error 
was  reduced  approximately  70S(  by  the  addition  of  the  compensator.  The 
design  procedure  was  based  on  the  Crossover  Model  for  the  human  pilot 
as  developed  by  Systems  Technology  Incorporated  of  Hawthorne,  California. 
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A DIFFERENTIAL  THRUGT  CONTROLLER 


FOR 

AIR  CUSHION  LANDING  SYSTEM  AIRCRAFT 


I.  Introduction 


Background 

Aircraft  dependence  upon  the  wheel  for  ground  operation  has  often 
presented  great  problemo,  especially  with  heavy  aircraft.  The  surface 
used  for  taking>off,  landing,  and  taxiing  has  to  be  strong  enough  to 
accept  the  large  force  exerted  by  the  small  area  of  tire  contact.  This 
causes  considerable  expense  in  airfield  construction  and  reduces  the 
mobility  of  air  forces.  Also,  the  wheels  have  to  be  attached  to  the 
aircraft  with  strong  structures,  adding  considerable  weight. 

A solution  to  this  problem  has  been  Bought  for  the  past  five  or 
six  years,  in  the  form  of  the  Air  Cushion  Landing  System  (ACLS).  This 
apparatus  is  a replacement  for  the  conventional  landing  gear.  The  main 
piece  of  equipment  is  an  oval  tube,  called  the  trunk,  rrhich  is  attached 
to  the  bottom  of  the  aircraft  fuselage.  An  auxiliary  engine  on  the 
aircraft  fills  the  trunk  with  air  and  forces  the  air  out  through  holes, 
located  on  the  bottom  and  inside  of  the  trunk.  As  a result,  an  air 
cushion  is  created  between  the  trunk  and  the  surface  of  the  ground.  In 
this  way,  the  weight  of  the  aircraft  is  evenly  distributed  over  the  total 
area  of  the  air  cushion,  resulting  in  a lower  ground  contact  pressure 
than  is  possible  with  tires  in  the  same  amount  of  ground  area.  The 
aircraft  is  then  able  to  taxi  over  very  soft  surfaces,  such  as  dirt 
and  mud.  The  idea  for  this  system  came  from  the  hovercraft,  a water 
vessel  designed  in  the  1930's  in  England. 
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The  Bell  Aerospace  Company  of  Buffalo,  New  York  hat  designed  and 
flown  an  ACLS  on  a single-engine  aircraft.  Tests  with  this  system 
revealed  that  it  was  possible  to  go  over  mud,  ditches,  small  tree 
stumps,  and  other  objects  that  would  have  stopped  a wheeled  vehicle. 

Also,  taking-off  arid  landing  on  unimproved  strips  wan  possible.  The 
successful  UBe  of  an  ACLS  on  a small  aircraft  has  paved  the  way  for 
further  development. 

The  Air  Force  Flight  Dynamics  Laboratory  at  Wright-Patterson 
Air  Force  Base  is  presently  involved  in  the  testing  of  a twin-engine 
aircraft,  equipped  with  an  ACLS.  This  aircraft  is  a Canadian  "Buffalo" 
and  is  designated  the  XC-8A  by  the  U.S.  Air  Force.  Through  scale  model 
testing,  and  computer  runs  with  a mathematical  model,  the  design  has 
been  accented  for  flight,  and  the  aircraft  1b  presently  undergoing  an 
extensive  flight  testing  program.  The  XC-8A  is  an  important  step  in 
the  development  of  the  ACLS  for  use  on  other  and  larger  military 
aircraft. 

Problem 

A problem  that  has  been  discovered  in  the  testing  is  the  suscep- 
tibility of  the  aircraft  to  movement  in  crosswindB  while  taxiing. 

Since  there  1b  no  tire  contact  to  give  directional  stability,  the  air- 
craft can  easily  be  blown  about  by  light  to  moderate  winds.  The  initial 
solution  attempted  was  to  install  exhaust  ports  on  the  outer  sides  of 
the  trunk.  These  exhaust  holes,  also  called  puff  ports,  can  be  individ- 
ually controlled  by  the  pilot,  to  push  the  aircraft  in  one  way  or 
another.  However,  the  forces  from  the  puff  ports  are  only  adequate  in 
a very  light  wind. 
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As  a result  of  the  directional  control  problem,  attention  war; 
turned  to  a powerful  source  of  yaw  control,  differential  thrust.  The 
potential  difference  in  thrust  between  the  two  engines  is  the  greatest 
source  of  low  speed  yaw  moment  in  the  XC-8A  (Ref.  5:150).  The  engine 
thrust  can  be  changed,  or  even  reversed,  by  changing  the  pitch  of  the 
propellers  and/or  the  throttle  setting.  Since  many  aircraft  are  steered 
on  the  ground  by  movement  cf  the  rudder  pedals,  it  would  bo  natural 
to  control  the  differential  thrust  by  movement  of  the  rudder  pedals. 
However,  the  time-lags  associated  with  the  engines  and  the  pilot,  and 
the  large  inertia  of  the  taxiing  aircraft,  dictate  the  need  for  a 
compensation  system  in  the  differential  thruBt  controller. 

Objectives 

It  wan  desired  to  design  a control  Bystem  that  would  improve  the 
airc.  '.ft  heading  response  to  the  pilot's  input,  with  a random  heading 
disturbance  due  to  the  wind.  It  was  alBo  desired  to  tost  the  system 
with  and  without  the  controller  to  determine  the  improvement  obtained. 

Approach 

Developing  an  adequate  model  for  the  engines  and  the  aircraft 
was  the  firet  Important  part  of  the  project.  Then  the  design  portion 
was  accomplished  uBing  classical  control  theory  and  the  well  known 
Crossover  Model  for  the  pilot.  It  was  assumed  that  future  application 
of  the  ACLS  will  be  on  Jet-powered  aircraft,  where  the  engine  time-lags 
are  higher  than  on  the  conventional-powered  XC-8A.  Therefore,  larger 
lags  than  the  one  associated  with  the  XC-8A  engines  were  considered 
for  further  design  purposes.  The  simulations  of  the  controlled  and 
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original  system  were  performed  on  a digital  computer,  using  the  MIMIC 
simulation  language,  which  is  a means  of  performing  analog  type  sim- 
ulations by  digital  techniques.  MIMIC  was  used  because  the  program 
was  easier  to  set  up  than  an  analog  program,  and  the  time-lags  were 
more  accurately  simulated  on  a digital  computer  than  on  an  analog, 
where  a Pade  approximation  would  hav'j  to  be  used.  Finally,  a mean 
square  error  analysis  of  the  system  was  perfcr.  1. 
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II.  The  Plant  Model 

As  a preliminary  to  modeling  the  aircraft  and  its  engines,  it 
was  necessary  to  determine  how  the  pilot  interacted  with  the  aircraft 
in  the  taxi  problem.  Much  thought  and  research  wan  dedicated  to  this 
fundamental  tank  and  the  resulting  block  diagram  can  be  seen  in  Fig.  1. 

While  the  pilot  is  taxiing  the  aircraft,  he  is  involved  in  a 
multiloop  task.  There  are  two  loops,  one  in  which  the  pilot  is  con- 
trolling the  heading  of  the  aircraft  (¥),  which  is  referenced  to  the 
direction  of  the  taxiway,  and  a Becond,  outer  loop  in  which  he  is 
controlling  the  aircraft  distance  from  the  taxiway  centerline  (y). 

Each  loop  has  a random  input  disturbance  reused  by  the  wind.  The 
theme  for  this  raultiloop  approach  came  from  two  sources,  References 
19  and  21,  in  which  the  pilot  controls  the  pitch  of  an  aircraft  and 
the  distance  off  the  glide  slope  while  on  an  instrument  approach. 

An  analyaia  of  the  complete  multiloop  task  was  beyond  the  Bcope 
of  this  thesis i however,  it  was  believed  that  when  control  of  heading 
in  the  inner  loop  is  optimised,  the  outer  loop  response  will  be  satis- 
factory. Therefore,  the  inner  loop  vaB  investigated  and  a controller 
designed  to  improve  its  characteristics. 

The  block  diagram  of  the  heading  control  loop  can  be  seen  in  Fig. 
2,  where  the  task  appears  aB  a regulator  control  problem  and  the  com- 
manded heading  (H1  ) will  be  sero.  This  block  diagram  is  topologically 
c 

equivalent  to  the  block  diagram  shown  in  Fig.  3.  If  the  portion  of 
the  block  diagram  incide  the  broken  line  in  Fig.  3 is  considered,  the 
problem  becomes  a compensatory  tracking  task  and  can  be  analysed  by 
the  Crossover  Model  theory  presented  in  Chapter  IV. 
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The  first  part  of  the  plant  to  be  Btudied  was  the  engines.  The 
engines  receive  their  input,  called  the  power  lever  angle  (PLA),  from 
the  power  levers,  which  control  thrust.  In  the  lower  thruQt  ranges, 
used  while  taxiing,  the  thrust  is  changed  by  a combination  of  changing 
both  power  and  the  blade  angle  of  the  propellers , The  blade  angle  of 
the  propellers  is  changed  by  the  beta  control  Bystem,  which  is  a mod- 
ification of  the  original  engine  control.  The  engine  output  is  thrust, 
which  converts  directly  to  yawing  moment  (N^)  by  multiplying  thrust 
times  the  distance  of  the  engine  from  the  aircraft  center.  Therefore, 
the  transfer  function  needed  for  the  engines  is  N^/PLA. 

Ah  mentioned,  the  engines  on  the  XC-8A  have  a beta  (propeller 
blade  incidence  angle)  control  ttyutem.  This  in  a modification  of  the 
original  engine  control  and  was  incorporated  because  the  original 
control  gave  only  a gross,  stopped  variation  of  thrust,  at  the  lower 

H 

power  setting.  The  installation  of  the  ACLS  o.u  the  aircraft  dictated 
the  need  for  a vernier  adjustment  of  thrust  while  taxiing.  One  result 
of  adding  the  bets  control  system  was  the  reduction  of  the  pure  time- 
delay  of  the  engine  response  to  about  0.2  second.  This  is  considered 
to  bo  a very  .low  delay  time,  as  the  delay  times  for  turbojet  and  turbo- 
fan engines  ere  on  the  order  of  one  second.  (This  time-lag  information 
was  obtained  from  discussions  with  Dr.  George  Kurylowieh  and  with 
Mr.  Elisha  Raohovitaky  of  the  Air  Force  Flight  Dynamlos  Laboratory, 
Wright-Patterson  Air  Force  Bass).  Since  larger  delay  times  may  be 
encountered  in  future  ACL8  airaraft,  it  was  decided  to  perform  the 
design  study  for  values  of  engine  reeponee  pure  time-delay  of  0.2, 
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0 . U , 0.6,  0.8,  and  1.0  seconds.  In  the  ground-based  simulation  of  the 
XC-8A  with  Air  Force  test  pilots,  the  differential  thrust  control  vsg 
simulated  by  a direct  mechanical  hookup  from  the  rudder  pedals  to  the 
power  levers.  The  engine  reuponoe  time-delay  vna  varied,  and  when  it 
was  increased  to  ona  second,  the  pilots  could  not  effectively  control 
the  simulator  while  taxiing  (Ref.  10). 

The  remainder  of  the  engine  dynamics  were  simulated,  in  this 
report,  as  a first  order  lag.  From  the  thrust-versus-time  response 
curve  of  Fig.  h , it  appears  that  the  thrust  response  to  a step  input 
of  PLA  is  velj.  approximated  by  a first  order  lag  with  a time  constant 
of  0.9  second.  Both  the  thrust  response  and  the  laeal  first  order 
lag  curves  are  drawn  for  comparison.  The  maximum  steady-state  thruat 
of  2000  pounds  is  used  because  this  ii  the  limit  of  the  beta  oontrol 
system  and  is  considered  sufficient  for  taxi  purposes.  The  maximum 
amount  of  yawing  moment  that  can  be  obtained  from  i 2000  pounds  of 
thrust  is  69,000  foot-pounds  (Ref.  5:191).  Therefore,  the  complete 
engine  transfer  function  is 

PLA  1 + 0.98  (1) 

where  Ng,  is  in  foot-pounds,  PLA  la  in  degrees,  and  1 a is  the  engine 
response  pure  time-delay. 

The  aircraft  model  required  the  use  of  many  engineering  approxi- 
mations. The  effect  of  the  vertical  stabiliser  on  the  yaw  moment, 
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due  to  the  and  stability  derivatives,  was  the  only  aerodynamic 
R $ 

effect  considered.  This  assumption  is  Justified  by  the  low  taxi  speed. 
The  aerodynamic  effect  considered  was  necessary  because  it  accounts 
for  the  heading  input  disturbance.  There  is  one  aerodynamic  effect 
neglected  which  may  have  some  significance,  the  rolling  moment  caused 
by  propeller  wash  over  the  wings.  It  was  assumed  that  this  can  be 
eliminated,  if  undesirable,  by  means  of  automatic  spoiler  or  aileron 
control.  It  was  also  assumed  that  the  pilot  keeps  the  forward  speed 
of  the  aircraft  constant  at  twenty  knots  ( 33.76  feet/second)  by  con- 
trcllin^  power  and  brakes.  (The  brakes  are  expandable  pillows  on  the 
bottom  of  the  trunk,  with  attached  Bkid  pads).  Finally,  the  yaw 
moments  due  to  the  brakes  and  the  puff  ports  were  neglected. 

Wit.i  the  necessary  assumptions  made,  the  model  for  the  aircraft 
can  be  derived.  The  model  must  account  for  two  inputs  and  two  outputs, 
so  each  Input/output  pair  will  be  discussed  separately.  The  engines 
deliver  a yaw  moment  input  to  the  aircraft,  which  1b  converted  into 
aircraft  heading.  The  model  for  this  conversion  process  is  a simple 
inertia  transfer  funotiont 


zz 


(2) 


where  H*  is  the  heading  due  to  the  engines  and  I is  the  moment  of 

£4  ZZ 

inertia  about  the  aircraft  z-axis. 
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The  other  input  to  the  aircraft  is  the  wind  disturbance,  and  the 
aerodynamic  effect  of  the  wind  on  the  heading  of  the  aircraf.  is 
mainly  due  to  the  vortical  stabilizer.  The  geometry  of  the  wind 
effect  problem  can  be  seen  in  Fig.  5.  The  following  equation  accounts 
for  the  yaw  due  to  wind: 


N 


W 


h ov2sb 


( V * c» 


R 


(3) 


where  p iB  the  atmospheric  density,  S 1b  the  wing  area,  b is  the  wing- 
span, and  R is  the  yaw  rate.  The  relative  wind  velocity  is  found  from 
the  relationship 


V - 


(>0 


where  and  are  the  relative  wind  velocities  along  the  X and  Y 
aircraft  axeB,  and  the  sideslip  angle  1b  defined  as 


e - Y 
V 


Substituting  equations  (U)  and  (5)  into  (3)  yields 

»» ■ | sevY  * °»Bb  i 


and  since 


NW  - Vv 


(5) 


(6) 


(7) 
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and 


R - ¥ (8) 


(assuming  zero  roll  and  pitch)!  then 


Wv  • K.#Tl/M  * Wvx*yriv  (9) 


where  Kj  ■ *s  pSbCjj  and  K2  ■ pSb2Cjj 
0 R 


From  the  geometry  of  the  problem  in  Fig.  5* 


Vx  - 33.78  cos  Yw  + Vg  Bin 

VY  ■ -33.T8  Bln  *„  ♦ Vg  co.  V,  (10, 

Assuming  that  the  heading  angle  will  be  kept  Bmall  by  the  pilot , the 
small  angle  approximations  can  be  applied  and  (10)  becomes 

Vx  « 33.78  + Vgfw 

VY  - -33.78  ¥w  + Vg  (11) 

It  is  desired  to  linearize  equation  (9)  by  a Taylor  Series  expansion. 

• •• 

The  operating  point  is  defined  as  V ■¥  ■¥  ■¥ ■ 0 and,  us  a result 

go  wo  wo  wo 

of  (11),  Vxo»  33.78  and  VY0-  0.  When  equation  (ll)  is  substituted 
into  (9),  the  result  is  of  the  form 


■ f<wv  ■ «<vw 


(12) 
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The  Taylor  Series  expansion  is 
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and,  by  the  Chain  Rule, 


(V„  - V ) + H.O.T, 
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At  the  operating  point: 


4-  ■33.T8K,.|t  -0.  f£  ■ 33.78  K,, 


av  av 

S=r  - -33-78,  and  ■==*■ 

3V  3Vg  34, 


1.  so  • 33.78  K2,  |f  » ~(33.70)2K1, 


and  !§■  ■ 33.78  K. . Ignoring  the  higher  order  terms  (H.O.T. ),  the 

avg 

Taylor  Series  expansion  about  the  operating  point  becomes 


!zs4v  - 33.78  K24w  - (33.78)2Kj4v  + 33.78  KjV 


(17) 
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When  equation  (17)  is  put  into  v.rannfcr  function  form  and  the  parameter 
values  from  Fig.  6 are  substituted,  the  result  is 

. 0.0007518 

Vg  s2  + 0.0756s  + 0.0254  (18) 

The  aircraft  acta  as  a second  order  filler,  with  & damping  factor  of 
0.24  and  a natural  frequency  of  O.lfi  radian  per  second.  These  numbers 
agree  with  the  damping  factor  and  natural  frequency  of  the  Dutch  Roll 
mode  for  the  same  stability  derivatives,  which  serves  as  a good  check 
on  the  mathematics.  The  block  diagram  of  the  plant  model,  with  numer- 
ical values  inserted,  can  be  seen  in  Fig.  7. 
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Engines: 

Propellers i 
Wing  Area: 
Wing  Span: 
c : 

Weights : 


2 x GE/T61*-ll* 

2 x UCAL  ST  6-73  (for  ACLS) 

Ham.  Std.  63-E60-15,  3 bladed,  H*.5  ft.  dian. 
9^5  sq.  ft. 

96  ft. 

10.3  ft. 

Max.  Takeoff  1*1,000  lb.  (structural  limit) 
Max,  Landing  39,100  lb.  (structural  limit) 

Low  Speed  Stability  Derivatives  (Ref.  12) 


0.22 


Cv  ■ 0.105 

Ns 


Inertia  (1*1.000  pound  A/C.  1*1.5#  c) 

I - 50861*2.0  slugs  ft2 

22 

Sea  Level  Air  Density 
p - 0.00237690  slugs/ ft ® 


Fig.  6 Relevant  Data  for  the  XC-8A  Model 


Fig.  7 Block  Diagram  of  the  Plant  Model 
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III.  The  Wind  Model 

Thu  wind  model  accounted  for  the  only  disturbance  to  the  plant 
and  vas  used  to  determine  the  effectiveness  of  the  controller.  Many 
different  wind  models  were  found  in  a literature  search,  among  them 
the  well-known  vonKarman  and  Dryden  wind  models,  but  they  were  all 
functions  of  the  aircraft  speed.  These  modelB  all  assume  the  aircraft 
to  be  flying  through  a gust  field  which  1b  frozen  in  time  and  do  not 
apply  very  well  to  a taxiing  aircraft,  where  the  gust  effect  ia  from 
a combination  of  temporal  and  spatial  fluctuations,  Since  a model  for 
t'hiB  was  not  available,  and  there  was  not  sufficient  time  to  collect 
data  and  analyze  it,  the  wind  was  approximated  as  a band  limited  white 
noise,  with  a bandwidth  of  one  radian  per  second.  Band  limited  white 
noise  iB  defined  to  be  a noise  with  a constant  power  Bpectral  density 
(PSD)  over  a finite  range,  as  seen  in  Fig.  8.  The  exact  cutoff  fre- 
quency was  not  critical  as  long  as  it  was  high  enough,  since  the  noise 
wsb  cut  off  at  a much  lower  frequency  aB  a result  of  filtering  by  the 
aircraft.  This  lower  cutoff  frequency  was  baaed  on  the  assumption 
that  no  structural  vibrations  were  to  be  modeled. 

The  wind  was  assumed  to  be  blowing  across  the  taxlvay  and  con- 
sisted of  the  gusting  effect  only.  The  steady  component  of  crosswind 
was  not  considered  important  to  the  problem  as  the  pilot  could  sense 
it  and  apply  a steady  correction.  However,  it  does  affect  the  dynamics 
of  the  aircraft  response  by  changing  the  coefficients  in  equations  (13) 
and  (l8).  As  the  steady  component  of  wind  increases,  the  guBt  effect 
on  the  aircraft  decreases  and  the  aircraft  response  becomes  overdamped. 
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Since  thia  ia  an  improvement  over  the  lightly  damped  case*  it  was 
decided  to  eliminate  the  steady  component  of  croaawind  and  deBign  for 
the  worBt  case.  Thus,  a mean  of  zero  was  selected  for  the  wind  gusts. 
The  deteimlnation  of  the  variance  presented  a more  difficult  problem. 
Since  the  aircraft  was  restricted  from  flying  when  the  wind  gusts 
exceeded  a certain  level,  the  intensity  of  the  maximum  guata  encountered 
was  under  human  control.  Assuming  that  the  average  guBta  encountered 
would  be  ten  knots,  or  ± five  knots  about  the  mean,  a standard  devia- 
tion of  five  knots  (8.U5  feet  per  second)  vaB  selected  for  the  wind, 
and  the  resulting  power  spectral  density  can  be  seen  in  Fig.  8.  The 
calculations  for  this  figure  can  be  found  in  Appendix  D. 
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IV.  The  Pilot  Model 

T'ick.cround 

In  order  to  analyze  the  heading  control  problem  and  then  to  design  ! 

a control  system  to  correct  it,  It  war  necessary  to  have  a model  for 
the  pilot  in  the  inner  control  loop  of  Fig.  1.  Extensive  documentation 
was  found  on  the  subject  of  pilot  modeling,  especially  in  references  1 

13,  lit,  15,  and  16;  and  the  bibliography  in  reference  16  provided  a 
good  source  for  a literature  search.  The  model  selected  was  the  moat 
widely  accepted  one  in  the  field  of  control  systems  engineering,  the 
Crossover  Model.  Since  no  deterministic  method  waB  found  to  write 
down  a mathematical  model  for  the  human  pilot,  many  engineering  Judgment* 
had  to  be  made  in  the  selection  process.  The  deBign  and  final  analysis 
of  the  design  were  based  on  these  Judgments,  which  could  not  be  vali- 
dated by  experiments  with  human  pilot's,  due  to  lack  of  time  and 
equipment . 

The  Crossover  Model  was  derived  from  a quasi-1 inear  representation 

i 

of  the  non-.linear,  time-varying  human  pilot.  The  pilot  can  be  repre-  j 

f 

sented  over  a wide  range  of  frequencies  by  a transfer  function  plus  S 

a remnant.  One  form  of  the  transfer  function  is  (Ref.  16) 


TlJui+  l\ 

TjJu)  + 1 J 


1 \ 

\T>+1/ 


(19) 


where 


e-j0JT 


» gain 

■ fixed  time-delay  due  to  conduction  time  of 
various  subsystem  elements 
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+ 1 m characteristic;  and 

T^Ju.'  + 1 values  depend  on  form  ol’  controllt=d  element 


TkJu>  + 1 
T^Jui  + 1 


“ terms  which  arc  uned  to  describe  very  low 
frequency  phase  data 


and  the  term  in  brackets  accounts  for  neuromuscular  system  elements. 
Thin  is  called  the  preciaion  model  and  accounts  for  the  major  sub- 
systems in  the  behavior  of  the  pilot.  Tho  variables  can  change  from 
one  person  to  the  next  and  can  even  change  for  tho  uume  person  at  a 
single  task  as  futigue,  boredom,  motivation,  and  other  human  ta^tors 
change. 

The  remnant  is  a term  used  to  account  for  the  difference  between 
the  output  of  the  human  pilot  and  the  output  of  the  model.  It  ie 
represented  by  a power  spectral  dennity  of  a noise  injected  into  the 
system  at  the  pilot's  output,  and  is  mostly  caused  by  non-linear, 
time-varying  pilot  bohavior.  The  remnant  term  has  been  refined  over 
the  years  but  la  still  difficult  to  model.  In  this  project,  the 
remnant  was  not  considered  important  at  the  frequencies  of  interest 
(about  one  to  five  radians  per  second). 


The  Crossover  Model 

The  Crossover  Model,  as  discussed  here,  applies  to  the  compensatory 
tracking  task,  as  seen  in  block  diagram  form  in  Fig.  7.  In  the  com- 
pensatory task,  the  pilot  rsceives  a visual  signal,  equlvolsnt  to  tha 
difference  between  the  forcing  function  and  the  output  of  the  controlled 
element,  and  attempts  to  minimize  this  error  signal.  For  the  Crossover 
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Fig.  9 Tho  Compensatory  Tracking  Task 


Modal  to  be  valid,  the  forcing  function  nuat  be  a random  appearing 
Input  and  the  signal  to  the  pilot  must-  bo  visual.  The  effeota  of  other 
pilot  cues,  such  as  motion,  have  been  studied  elsewhere  (Ref.  lB). 

By  experimental  observation,  the  characteristic#  of  the  open  loop 
transfer  function  of  the  compensatory  system  have  been  found  to  be 
elmllar  l'or  many  different  controlled  elements.  The  open  loop  log- 
mngnltudo-vereua-log-fraquency  (Bode)  diagrams  all  have  a straight 
line  with  a slope  of  minus  twenty  decibels  (db)  per  decade  in  the 
region  where  the  log  magnitude  crosses  the  aero  db  point.  This  is 
where  the  term  "Crossover  Model"  origtnutee.  In  othor  words,  in  the 
crossover  region  of  frequency,  the  open  loop  transfer  function  is 

Vc  " "^s  (SO) 

whero  is  the  crossover  frequency  and  t is  the  time  delay  due  to 
both  the  pilot  and  the  controlled  element.  It  is  obvious  that  ths 
Crossover  Model  of  the  pilot  is  dfpendant  on  the  form  of  the  con- 
trolled  element. 
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When  tha  controlled  element  tranofer  function  la  known,  the 
nodal  for*  tha  pilot  in  the  oroBuovcr  region  can  bo  found  from 

a.  cTT«“ 

Y . ~.JL_ 

K «yc  (ai) 

whera  t ia  tha  time-delay  associated  with  tha  pilot  and  any  tlrau- 
delayn  in  the  controlled  element  are  not  oonaidercd.  Thla  ia  easily 
dona  mathematically,  but  tha  pilot  la  not  capable  of  efficiently 
becoming  a second  order  derivative  or  higher.  When  he  ia  saked  to  do 
ao  i hi  a opinion  of  the  controlled  element  neriuualy  deteriorates. 

Thla  opinion  can  be  objectively  rated  on  a aoale  from  one  to  ton, 
called  the  Cooper-Harper  Beale » which  la  well  known  to  teat  pilots 
and  aircraft  dettlgu  engineers  as  a ban  In  for  rating  the  handling 
Huai  itleit  of  al  re  raft  (Fig.  10). 

The  pilot'd  time  delay  la  dependent  on  the  controlled  element 
and  on  the  forcing  function  bandwidth.  Thla  time  delay  is  difficult 
to  evaluate  deterministically  as  it  can  be  improved  with  training 
and  varies  with  different  people.  However,  it  can  be  predicted  with 
reasonable  aoouraoy  by  (Fief.  13) 

tYc>  - (sc) 

whore  TfJ  depends  upon  the  controlled  element  and  At^  depend*  upon  the 
forcing  function  bandwidth,  The  term  At#  ia  approximately  0.00  uij . 
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Demands  on  the  Pilot 

Aircraft  C 
Characteristics 

Pilot  compennation  not  a factor 
for  desired  performance 

Excellent 
Highly  desirable 

Pilot  compennation  not  a factor 
for  deuirod  performance 

Good 

Negligible  deficiencies 

Minimal  pilot  compensation 
required  for  desired  performance 

Pair  - Some  mildly 
unpleasant  deficiencies 

Desired  performance  requires 
moderate  pilot  compensation 

Minor  but  annoying 
deficiencies 

Adequate  performance  requires 
considerable  pilot  compensation 

Moderately  objectionable 
deficiencies 

Adequate  performance  requires 
extensive  pilot  compensation 

Very  objectionable  but 
tolerable  deficiencies 

Adequate  performance  not 
attainable  with  maximum 
tolerable  pilot  compensation. 
Controllability  in  question 

Major  deficiencies 

Considerable  pilot  compensation 
lu  required  for  control 

Major  deficiencies 

intense  pilot  compensation  is 
required  to  retain  control 

Major  deficiencies 

Control  will  be  lost  during 
some  portion  of  required 
operation 

Major  deficiencies 

Pig.  10  Cooper-Harper  Handling  Qualities  Rating  Scale 

The  crossover  frequency  Is  another  subjective  variable.  It  is 

estimated  as  (Ref.  13} 

% " uo0(V  * ^o(“i)  ( 
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To  a first  approximation,  Auj^  iB  zero  for  the  forma  of  controlled 

elements  found  in  aircraft,  and  u can  be  found  from  the  condition 

co 

for  a neutrally  stable  system  at  the  crossover  frequency,  with  no 
input,  The  equation  is 


TU 


co 


V 

8 


(2M 


where  t is  the  sum  of  tq  and  the  controlled  elemant  time-lag,  if 
there  ie  one. 


The  ACLB  Pilot  Model 

As  seen  in  Fig.  5»  the  controlled  element  for  the  ACLS  aircraft 

ii 

Y „ 0.0043  e~"T» 

C a2  (a  + Ull)  (25) 

This  la  a very  difficult  element  for  the  pilot  to  control.  He  haa  to 
generate  a transfer  function  of  the  form 


Yp  « Kp  e"IT®  a(a  + 1.111)  (26) 

in  the  crossover  region.  This  means  that  the  acceleration  channel  of 
the  pilot  is  activated,  causing  his  zero  input  time-delay  (tq)  to  be 
about  0.8  second  (Ref.  l6i26)  and  the  Cooper-Harper  rating  to  be 
extremely  high  (about  nine  or  ten,  Ref.  16 i 37 ) • The  pilot  rating  alone 
is  sufficient  reason  to  reject  the  controlled  element  aa  undesirable. 
Since  it  la  only  a subjective  measurement,  the  system  was  analyzed 
from  the  Crossover  Model  theory. 
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The  lowest  value  of  the  controlled  element  time-delay  was  0.2 
□econd.  Adding  this  to  the  tq  value  of  0.8  Becond  and  using  equation 
(2l»),  the  crossover  frequency  was  found  to  be  1.57  radians  per  second. 
The  crossover  frequencies  for  the  other  controlled  element  time-delays 
of  0.4,  0.6,  0.8,  and  1.0  seconds  were  1.31,  1.12,  0.98,  and  0.8J 
radians  per  second,  respectively.  According  to  reference  one,  a 
realistic  value  of  crossover  frequency  for  handling  qualities  consid- 
erations la  about  2.0  ± 0.5  radians  per  second.  Thus,  the  crossover 
frequency  of  the  syatem  is  lower  than  desired,  Also,  the  lower  cross- 
over frequency  causes  a higher  syatem  rms  error  aB  the  error  is 
inversely  proportional  to  the  oquare  of  the  crorsover  frequency 
(Ref.  I6il5). 
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V . Control  Syotem  Design 


C 


Design  Philosophy 

The  control  system  design  was  based  on  the  idea  that  the  form  of 
the  controlled  element  determines  the  pilot 'b  describing  function. 

Also,  it  was  desired  to  reduce  the  engine  time-delay  as  much  as  possible. 
A general  formula  for  the  pilot's  describing  function  in  the  crossover 
region  is 


Y 


P 


Kp  a 


■ST, 


T.s  + 1 

Tji  + 1 


(27) 


where  and  Tj  are  the  amounts  of  lsad  and  log  that  the  pilot  must 
generate  to  compensate  for  the  controlled  element  dynamics.  The  lead 
(T^)  is  the  quantity  which  has  the  most  effect  on  the  pilot's  opinion 
of  the  system.  A reasonable  maximum  for  is  five  seconds  and  a 
graph  showing  the  effeot  of  lead  generation  on  the  pilot  opinion 
rating  is  seen  in  Fig.  11.  The  pilot  rating  scale  referred  to  is  the 
one  shown  in  Fig.  10, 


Fig.  11  Hating  Decrement  Due  to  Pilot  Lead  us  Inferred  from 

Handling  Qualities  Tests  (Ref.  1 U i SL  7 ) 

(Open  curve  due  to  uncertainty  of  T. ) 

h 
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According  to  Ref.  1,  the  boot  handling  qualities  are  obtained 

when  tho  value  of  1 is  Icon  than  about  1.5  seconds.  It  was  decided 
L 

to  make  the  value  of  for  this  design  study  equal  to  one  second.  As 
will  be  seen,  this  criterion  had  to  be  relaxed  somewhat  for  stability 
reasons.  For  the  pilot  lead  to  be  one  second,  the  controlled  element 
was  required  to  have  a factor  of  s + 1 in  the  denominator,  at  the  cross- 
over frequency.  Assuming  that  the  controlled  element  can  be  configured 
to 
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• («  + !) 
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then  the  pilot  describing  function  would  be 

Yj,  - Kp  o’“T«  ( s + 1 ) (29) 

in  the  crossover  region.  The  pilot's  lag  term  wan  not  presented  in 
aquation  (29)  because  it  was  assumed  that  the  Ideal  controlled  element 
would  have  no  lead  term  in  the  numerator,  This  lead  term  In  the  usual 
uauae  of  the  pilot's  adopting  a lag.  For  a controlled  element  of  thn 
form  of  equation  (2fl),  the  value  of  the  pilot's  pure  time-delay  with 
Kero  Input  (tQ)  is  approximately  0.33  seaond  (Ref,  13),  (Other  values 
of  time-delay,  ranging  from  0,16  second  to  0.50  second,  were  found  in 
the  literature.)  Hie  objective  of  the  control  system  design,  then, 
was  to  closely  approximate  Ky/sU  + i)  in  the  crossover  region. 

The  first  step  in  the  classical  control  theory  design  procedure 
was  to  change  the  controlled  element  into  a feedback  control  system. 
From  a study  of  thn  original  system,  it  was  determined  that  ths 
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quantities  available  for  feedback  were  the  heading  and  the  heading 
rate.  The  heading  rate  was  easier  to  detect  by  using  an  available 
piece  of  aircraft  hardware,  the  yaw  rate  gyro.  With  the  yaw  rate 
(heading  rate)  fed  back  to  the  input  of  the  engines,  the  second  step 
consisted  of  inserting  a cascade  compensator  into  the  plant  ahead  of 
the  engines.  The  resulting  block  diagram  for  the  controlled  system 
is  seen  in  Fig.  12. 


Design  Method 

The  block  diagram  of  Fig.  12  shows  a block  containing  a 1/s  (an 
integrator)  outside  of  the  control  loop.  Thin  Integrator  contributes 
the  s term  in  the  denominator  of  equation  (20),  the  desired  transfer 
function  of  the  controlled  element.  Thus,  the  objective  in  designing 
the  compensation  control  loop  (CCL)  was  to  make  it  asnums  the  closed 
loop  transfer  function 


CCL 


is. 

s ♦ 1 


(30) 


The  method  for  the  design  of  the  compensator  for  the  feedback 
aontrol  system  come  from  reference  22.  Tile  principle  was  to  equalise 
the  nystem  as  if  the  time-delay  waa  not  there  and  then  to  compensate 
for  the  time-delay.  The  time-delay  was  reduced  by  adding  lead  com- 
pensation to  the  system,  with  the  nereis  of  the  lead  anmpenttaior  out- 
side of  the  system  bandpass. 

The  first  step  taken  was  to  draw  the  hoot  Locus  of  the  CCL  with 
yaw  rate  feedback  and  no  oor.pensaUon  added.  In  Appendix  A,  an 
explanation  is  given  of  how  to  draw  the  Hoot  Locus  for  s feedback 
control  system  wi*h  time-delay.  The  hoot  Locus  plots  for  values  of 
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of  tHe  ttrmdTTic  Control  Loop  vith  Compensator  Added 


t equal  to  0.2,  0,k,  0.6,  0.8,  and  1.0  seconds  can  be  seen  In  Figures 
B.l  through  B.5  in  Appendix  B.  From  the  plots,  it  could  be  Been  that 
it  was  impossible  to  select  a gain  which  would  put  one  dominant  root 
at  -1.  If  one  root  was  put  at  -1,  then  another  root  would  be  more 
dominant  and  almost  at  the  origin. 

Therefore,  the  first  procedure  in  the  design  was  to  reconfigure 
the  Root  Locus  so  that  one  root  could  be  selected  at  -1  and  the  other 
.roots  would  be  less  dominant.  This  was  done  with  a compensator  of  the 
fora 


(31) 


where  was  the  amplification  factor  and  the  pole  was  ohosen  to  be 
ten  times  the  aero  for  mathematical  simplicity.  The  Root  Locus  plots 
for  the  system  with  this  lead  compensator  added,  for  the  various  values 
of  are  drawn  In  Figures  8,6  through  8.10  In  Appendix  n.  From  each 
of  then#  plots,  It  waa  possible  to  select  one  root  at  -.1. , with  the 
other  roots  being  loss  dominant.  For  the  ■ 1.0  second  case,  the 
second  root  was  very  close  to  the  first  at  -1,  but  this  was  corrected 
in  the  second  design  step, 

The  second  step,  designing  a multi-order  lead  compensator  to  re- 
duce the  pure  I. line-delay  of  the  engine  response,  was  more  difficult 
than  the  first.  The  transfer  function  for  the  compensator,  called  the 
time  delay  compensator  (TOC)  was 


31 


where  was  the  gain  constant,  a waa  the  order  (number  of  poles 
and  zeros),  and  2 and  p were  the  approximate  locations  of  the  zeros 
and  poles,  respectively.  Each  pole  was  set  equal  to  ten  times  the  zero 
to  keep  its  effect  negligible.  A pole  was  uned  with  each  zero  In  order 
to  keep  the  system  realizable.  The  phase  lead  contribution  of  the 
compensator  wan  approximately  n tan”*  — and  thin  phase  angle  had 

ft 

to  compensate  for  the  phase  lag  of  the  engine  time-delay.  The  two 
were  set  equal  at  the  crossover  frequency  to  yield 


n tan 


■ W T 

0 a 


(33) 


with  the  restriction  that  the  tero(s)  should  lie  outside  of  the  cross- 
over frequency.  Reference  {??  suggests  a limit  on  the  number  (n)  of 
zeros  to  keep  the  amplitude  response  of  the  system  respectably  low  for 
higher  frequencies,  The  number  of  zeros  was  net  at.  t.wo  In  this  rase 
because  the  system  had  two  poles  near  the  lower  frequencies  or  Interest, 
lly  using  the  value  of  0.33  second  for  the  pilot's  zero  input,  time-delay 
and  assuming  that  the  engine  time-delay  would  bo  reduced  to  zero,  the 
crossover  frequency  found  from  equation  (fit)  was  14,76  radians  per  sec- 
ond. Using  this  value  for  the  crossover  frequency  In  equation  (33), 
and  by  setting  n equal  to  two,  the  approximate  locations  of  the  poles 
anil  zeros  were  found.  As  the  value  of  ■)  Increased,  it  was  assumed 
that  the  crossover  frequency  would  decrease  because  of  Inability  to 
eliminate  all  of  the  phase  lag  due  to  i^.  The  assumed  values  of  cross- 
over frequency,  order  of  the  compensator,  and  selected  locations  of 
the  pules  and  zerus  for  the  five  1 values  can  be  seen  In  Table  1.  As 
can  be  seen,  it  was  necessary  to  decrease  the  crossover  frequency  to  the 
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Time  Delay  Comoensator  Parameter  Values 


Parameter 


ml ni mum  locoptnblr  an  the  value  of  i Increased,  to  heap  r,  larger 

a 

than  ui0. 

With  the  entire  compenual  Ion  ay  stem  dnelgned,  except  for  the 
amplifier  gains,  the  next  step  was  to  draw  the  Hoot  hocus  plots  for 
each  sy u t piii  mu  that  the  tens  It  l v I ty  , or  open  loop  gain,  could  he 
selected.  The  hunt  litmus  plots  for  the  five  systems  are  drawn  in 
figures  h.ll  through  H.lt  to  Appendix  11.  from  these,  the  sensitivity 
l.liat  caused  the  posit  Ion  of  the  dominant  root  to  he  at  -I  was  unite  bud. 
Tills  sensitivity,  K,  Is  actually  the  product  K^K.pi^.to.nuli  101) , Using 
this  sensitivity,  s closed  loop  Mode  plot  of  the  t'Oh  was  drawn.  This 
was  compared  to  I, he  Undo  plot  of  a tranufar  function  eipinl  to  TjS  + l/ 

Tj  a + I to  determine  (.he  amount  of  lead  and  lap  that  the  pilot  would 
Utue  t.o  generate  to  compensate  for  the  system,  finally,  the  total 
system  was  checked  for  stability  by  adding  the  freiiuehcy  response 


£ curves  of  the  pilot  modal,  the  CCL,  and  the  integrator.  In  two  oases, 

for  t equal  to  0.8  and  1.0,  it  van  necessary  to  lower  the  sensitivity 
of  the  CCL  to  make  the  total  system  stable.  Since  each  of  the  values 
of  engine  time  delay  required  a different  design  analysis,  they  will 
be  dijcusaed  separately. 

Tha  CCL  closed  loop  Bode  plot  for  equal  to  0.2  second  is  seen 
! in  Fig.  13.  The  gain  used  was  1018.29  because  that  gain  caused  the 

i 

) dominant  root  to  lie  at  -1.  The  log  magnitude  plot  closely  resemblea 
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a plot  of  the  transfer  function  1/s  + l,  which  is  drawn  as  a broken  line 
in  Fig.  13.  The  two  log  magnitude  curves  are  similar  up  to  u>  • *4.0. 
Fiinoo  the 'pilot  attempt!!  to  cancel  the  dynamics  of  the  CCL  in  th* 
crossover  region,  his  transfer  function  is 


Yp  - Kp  o_#T®  (b+1)  (3U) 

which  Is  exactly  as  desired.  The  pilot  gain  (Kp)  is  the  gain  needed 
to  raise  the  total  system  open  loop  log  magnitude  curve  so  that  it 
oroBsoa  the  Kero  db  line  at  the  crossover  frequency.  The  crossover 
frequency  was  found  from  the  condition  for  neutral  stability  with  *#ro 
Input.  At.  the  crunaover  frequency  of  u » 3.6,  the  pilot's  phene  angle 
(assuming  the  Kero  Input  pilot  tlnm-ti*  toy  wan  0.33  second)  va»  the 
t'Cii  phase  at ik it’  wen  ~9h”,  and  the  Integrator  phase  angle  waa  **90°.  The 
pilot's  gain  wan  net  equal  to  the  ornnso.-er  frequency  '11  this  and  the 
following  cases. 

The  01?!  1 c lotted  hu  p llodr*  plot  for  1 « 0,I|  second  Is  seen  111 

M 

Fig.  1 h . The  gain  used  was  HKjH.OT  and  till  a caused  the  dominant  foot, 
to  lie  at  -I*  The  log  magnitude  plot  clone  y resembles  a plot  of  the 
transfer  function  0,  S3  sH/e  + l , which  Is  drawn  as  a broken  line  In 
Fig . |li.  The  two  lug  magnitude  curves  are  similar  up  to  ui  * 1.0. 

The  pilot's  transfer  function  In  title  ona«  is 
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The  pilot  now  has  to  generate  a lag  of  0.  si  second  In  sddltlun  tc:  n 
lead  of  one  second.  The  crossover  frequency  Was  found  to  h*  J.i'3  radians 
per  second.  At  l.lilti  fi  equerr  y,  * lie  pilot's  phane  angle  was  the 

(H’l,  phene  angle  was  and  tre  integral  ir  phase  angle  was  -90". 
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Fig.  ,1  It  C ,1(18  0(1  IH'ni*  Plot,  of  001. , ■ ().!»  Boooml 

(firohtm  lim»  !«  otjul valwnt  transfer  function 
l on  inagni  tuda) 

Vli*  f.'fJi.  I udwil  luuii  Oniic  plot  for  i » O.fi  Booouil  In  mean  in 

Fig,  1!).  The  gain  u«#i1  in  thin  cm*  was  lWfi . 6 , can  a trig  the  dominant 

root  to  U«  at  >1.  'Hits  lug  magnitude  plot  roaomlileN  a plot  u(  th* 

traowfer  funotlon  i),67»U/a+l,  wlilnli  is  drawn  an  a broken  line  in 
« 

Fig.  lt>i  up  to  ui  ■ *'.l.  Ilia  pljot'a  transfer  funuliim  for  thia  caa«  in 
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The  pilot  now  ha*  i.o  generate  a lead  nf  one  aeoond  and  a lag  of  0.6T 
HoiMiinl.  The  Di’oBauver  frequently  and  pilot  gain  were  found  to  be  2,0 
rad  I ana  per  ueooinl,  At  this  frequency,  the  pi  lot 'a  phene  audio  waa 
-*.'0®,  the  (HU,  phase  angle  waa  -6un,  and  the  Integrator  phase  tingle 
waa  -90®, 

The  COL  u t used  loop  llutle  plot  for  i ■ 0.(1  srcnond  la  neon  In 
Pig,  Hi.  The  gain  weed  was  715. Mi  and  waa  the  gain  required  to  plaoe 
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Fig.  .1,6  Closed  Loop  Bods  Plot  of  CCL,  * 0.8  Leoond 

(Broken  lino  Is  equivalent  transfer  function 
log  magnitude) 

the  dominant  root  at  -0.6.  A higher  value  of  gain  aauned  the  system 
to  be  unstable  by  giving  It  a negative  gain  margin,  so  it  was  decided 
to  require  the  pilot  to  generate  a little  more  lead.  With  the  gain 
value  used,  the  gain  margin  vae  1.13  db  at  1.9  radians  per  seoond, 
with  the  pilot  t.lme-delay  Bet  at  0.9  second  as  explained  in  Chapter  VI. 
Hie  log  magnitude  plot  alosely  resembles  a plot  of  the  transfer  function 
e+l/1.67s+l,  which  is  drawn  as  a broken  line  in  the  figure,  up  to 
oi  ■ 1.8.  The  pilot's  transfer  funotion  is  now 
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(37) 


The  increase  in  the  pilot  lead  over  the  three  previous  cases  was 
undesirable,  but  should  not  Berioualy  degrade  the  pilot  rating,  and 
ia  still  a significant  improvement  over  the  uncompensated  syi  tem. 

The  crossover  frequency  whb  found  to  be  1.7  radians  per  second,  where 
the  pilot's  phase  angle  was  -20°,  the  CCL  phane  angle  wa  ■ -70°,  and 
the  integrator  phase  angle  was  -90°. 

The  CCL  closed  loop  Bode  plot  for  * 1.0  second  is  neen  in 
Fig.  17.  The  gain  used  was  1*76.7^,  because  higher  v&lues  of  gain  caused 
instability  in  the  system  by  causing  the  gain  margin  to  be  negative. 

With  the  value  of  gain  used  the  gain  margin  vbb  0.30  db  at  u>  ■ 1 . 7 
radians  per  second  and  the  dominant  root  was  positioned  at  -0.5.  The 
log  magnitude  plot  closely  resembles  a plot  of  the  transfer  function 

4 

0.9s+l/2a+l,  which  is  drawn  as  a broken  line  in  Fig.  17,  up  to  id  ■ 1.5. 
The  pilot's  transfer  function  is  now 
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(38) 


The  pilot  is  required  to  generate  a lead  of  two  seconds,  whloh  Is 
undesirable,  but  the  lead  cannot  be  reduced  if  the  system  is  to  remain 
stable.  The  amount  of  lead  required  should  not  be  a problem  to  the 
pilot  and  is  still  an  improvement  over  the  uncompensated  system.  The 
crossover  frequency  was  found  to  be  1,5  radians  per  second,  where  the 
pilot's  phase  ungle  was  -10°,  the  CCL  phase  angle  was  -80°,  and  the 
integrator  phase  angle  was  -90°. 
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Fig.  17  Closed  Loop  Bode  Plot  of  CCL,  t ■ 1.0  Second 

& 

(Broken  line  is  equivalent  transfer  function 
log  magnitude} 

A summary  of  the  data  discussed  in  all  five  cases  Is  presented 
in  Table  II.  As  the  engine  time-delay  was  increased,  the  Bystem 
performance  could  not  be  kept  as  high  as  desirable,  but  it  was  still 
a great  improvement  over  the  Bystem  without  the  compensator  added. 
The  improvement  was  mainly  because  the  pilot  had  only  to  generete  a 
small  amount  of  lead  rather  than  an  acceleration  form  of  dynamics. 
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The  simulation  van  ncaompl i shed  liy  uning  the  CDC  C’YBEtt/T^  digital 
computer  nt  the  Alt*  Forou  Institute*  of  Techno] ogy , Wright.-lVitterson 
AFB,  Ohio,  and  the  MIMIC!  programming  language  written  by  Mr.  ft.  .T, 
Hansom  and  Mr.  H.  ft.  Peterson  at  Wright-Pat.  ter  son  AFB,  Ohio.  A 
revision  to  the  original  program  wan  made  by  Prof.  C.  W,  Richard  of 
the  Mathematics  Department  at  Aft IT,  which  enables  one  to  obtain  plots 
of  the  computed  data.  MIMIC  is  a programming  technique  by  whioh 
analog  type  nlmulations  can  be  performed  on  the  digital,  computer, 
using  « fourth  order  Runge-Kutta  Integration  routine.  The  step  slse 
can  be  left  to  vary  ao  that  the  local  relative  error  Jo  kept  below 
5 x 10“6,  or  it  can  be  bounded  by  DTMAX  and  DTMIN.  Each  MIMIC  program 
allows  the  use  of  up  to  five  user  written  subroutine*,  written  in 
Fortran  language.  This  feature  wno  very  useful  as  the  time  delay 
function  in  the  MIMIC  language  used  an  excessive  amount  of  central 
processor  time,  whioh  caused  the  program  to  run  for  about  eight  times 
the  amount  of  time  that  the  system  would  bo  actually  operating.  A 
time  delay  subroutine  was  written  in  Fortran  language,  specifically 
for  thin  simulation,  whei ' the  integration  Btep  elite  was  fixed,  and 
the  program  ran  for  about  one  third  of  the  actual  system  operation 
tlmu . The  result  was  » considerable  savings  In  time  and  money. 

Two  separate  MIMIC  programs  ware  written,  one  for  the  uncompen- 
sated system,  and  one  for  the  compensated  system.  These  prograae 
can  bo  seen  In  Appendix  C and  an  explanation  of  the  programming 
technique  can  be  seen  in  Ref.  3. 
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Mind  Ouflt.  Simulation.  Thu  representation  for  the  wind  km 
previously  described  a*  a bond  limited  white  noise.  This  noise  was 
simulated  in  the  program  by  panel  tig  a train  of  pulses , whose  width 
was  0.2  second,  through  a firm!  order  filter.  The  amplitude*  t..f  the 
pulses  wnn  a distributed  random  variable,  obtained  from  the 

normal  random  number  generator  in  the  MIMIC  program.  'Hits  generator 
was  called  every  0.2  second,  causing  the  pulse  amplitude  to  change 
at  0.2  second  intervals.  The  method  of  representing  a band  limited 
white  noise  by  a pulse  train  is  described  in  Appendix  0.  The  inte- 
gration step  else  was  set  at  0.01  second  as  a compromise  between 
computation  time  and  accuracy  and,  therefore,  Shannon's  Sampling 
Theorem  (Ref,  9)  was  satisfied.  This  theorem  states  that  th«  sampling 
frequency  should  be  at  least  twice  the  pulse  frequency  in  an  informa- 
tion system.  The  pulses  were  passed  through  a first  order  filter  to 
nmk*  certain  that,  the  powsr  spectral  density  (POD)  had  a bandwidth 
of  one  radian  per  seoond.  If  S(in)  is  the  PSD  of  the  input  to  the 
filter  and  S(out)  in  the  output  POD,  then  according  to  Ref.  1"(, 


1 

1 

* 


S(out)  ■ | H(Jui)  | * fi(in)  (39)  | | 

II 

rv  r i unction  j.  • ?•«/ juitj.,  ?•*  attMiiuatos  cue  ruu  uy  turps  uu  mi.  w - i.w,  i 

i 

so  the  bandwidth  of  the  filter  is  one  radian  per  second. 

Pilot  Slmul atlon.  The  pilot  model  used  was  the  simple  Crossover 
Model  with  no  low  frequency  phase  data  terms  or  neuromuscular  terms,  ; 

un  explained  In  Chapter  VII.  The  crossover  frequency  and  the  resulting 
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pilot  fit  in  were  calculated  from  the  neutral  stability  criterion,  using 
the  aero  input  pilot  time  dolny,  According  to  the  .literature,  the 
phene  margin  llmt  occurs  in  the  use  of  (he  Croatiover  Model  results 
from  n ileereane  in  the  pilot  f liim-delay . The  only  analytical  meana 
('mind  to  mince  Uiin  wnn  hy  equation  (Pi’),  where  Lhr  time-delay 
decrease  wan  due  to  the  forcing  function  bandwidth.  In  this  case,  the 
bandwidth  of  the  disturbance  was  too  small  (about  O.J?  radian  per 
second)  to  have  any  effect  on  the  pilot  time-delay.  However,  the 
delay  could  also  be  reduced  by  such  intangibles  as  motivation,  train- 
ing, and  skill,  for  purposes  of  the  simulation,  it  was  assumed  that, 
for  the  uncompensated  system,  the  pilot  time-delay  would  be  reduced 
from  0,(1  second  to  0.f>  second,  and  for  the  compensated  system,  from 
0,  S3  second  to  0.5?  second.  Thane  value#  are  comparable  to  value*  of 
time -da lay  found  by  experiment  in  Inferences  13  and  16. 

Error  Analysis.  The  operation  of  the  system  wan  analysed  by 
'computing  the  mean  square  error  as  defined  in  Hnf.  SO.  tf  « is  the 
d! f fore  lies  between  the  output  and  the  input  in  a control  system, 
the  mean  aquar*  error  can  be  defined  a* 


e*  " -A™  / ** It  (40) 

j 0 

where  T is  the  elapsed  time.  This  computation  was  performed  as  part 
of  the  MIMI0  program.  Ten  runn  were  made  for  each  value  of  engine 
time-delay,  each  run  having  a different  random  noise  input,  and  the 
runs  were  made  for  both  the  uncompensated  and  the  compensated  systems. 
The  average  of  the  mean  square  error  was  calculated  for  each  sot  of 
ten  runs  so  that  a wide  statistical  spread  would  be  obtained. 
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Two  main  rcMUlte  were  obtained  from  the  simulation.  Plots  were 
drawn  of  the  Input  heading  dlnturbance  and  output  heading  vnrxus  time 
for  both  1 hr*  um'itni|u'hi«rttt>i1  mid  the  iinmpenrmtrd  system,  it  each  value 
of  engine  time-delay.  On  each  of  the  ten  plots,  shown  in  Figures  18 
through  27,  the  input  and  output  angles  are  scaled  on  the  y-axis  and 
time  is  scaled  on  the  x-axis.  The  corresponding  uncompensated  and 
compensated  system  response  plots  for  each  value  of  engine  time-delay 
show  that  Un*  pilot  modal  can  follow  the  Input  disturbance  more 
closely  with  the  compensator  added.  An  explanation  of  the  evident 
oscillation  of  trio  output  is  presented  in  Chapter  VII.  l’ha  plots 
were  nil  drawn  for  the  same  random  input  heading  disturbance. 

The  second  result  obtains!  was  the  collection  and  averaging  of 
the  mean  square  error  data,  Ten  runs,  each  with  a different  random 
input,  were  made  for  each  system  at  each  value  of  engine  time-delay, 
The  mean  square  error  at  .100  soconds  was  averaged  over  the  ten  rune 
for  each  case  so  that  the  uncompensated  and  compensated  systems  could 
be  compared.  The  tabulation  of  the  average  mean  square  errors  of  ths 
1.00  second  runs  is  seen  In  Table  III,  In  all  oases,  the  average 
mean  square  error  for  the  compensated  system  was  decreased  by  about 
'jfO#  of  the  uncompensated  system  value. 
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VII,  Conclusions  and  Recommendation* 

Conclusions 

Control  Systems  have  been  designed  to  improve  the  heading  response 
of  the  taxiing  ACLS  aircraft  to  pilot  input  of  differential  thrust  for 

five  different  value*  of  engine  time-delay,  Theoe  designs  were  baaed 
on  the  Crossover  Model  for  the  pilot  and  were  tested  by  using  the  name 
model,  so  they  are  dependent  on  the  model's  aocuracy.  However,  the 
Crossover  Model  has  been  in  use  for  almost  ten  years  and  nothing  was 
found  in  the  literature  to  dispute  it. 

In  the  graphs  showing  the  simulation  input  and  output  versus 
time,  some  high  frequency  oscillation  is  evident.  This  oscillation, 
it  is  believed,  would  be  eliminated  by  a human  pilot  due  to  the  high 
frequency  filtering  action  of  his  neuromuscular  reactions.  This 
filtering  is  modeled  by  the  neuromuscular  terms  in  equation  (19), 
which  were  not  used  in  this  project.  The  use  of  the  neuromuscular 
terms  requires  an  in  depth  study  of  physiology,  and  that  was  considered 
beyond  the  scope  of  this  thesis.  However,  according  to  Ref.  13,  the 
neuromuscular  frequencies  are  usually  in  the  region  of  from  10  to  15 
radians  per  second. 

The  average  mean  square  error  values  are  all  reduced  by  about 
70/S  with  the  addition  of  the  compensator.  The  average  mean  aquare 
error  for  the  compensated  system  with  on  engine  time-delay  of  0.8 
second  is  about  the  same  as  the  error  for  the  uncompensated  system 
with  a time-delay  of  0.2  second.  Similarly,  the  error  for  the  com- 
pensated system  with  a delay  of  one  second  compares  with  that  of  the 
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unoompensated  system  with  a delay  of  0.1*  ueeond.  Whether  the  compen- 
sated syatom  moan  square  errors  are  small  enough  can  only  be  determined 
by  experiments  with  human  pilots. 

Even  without  the  improvement  in  average  mean  square  error,  BUb- 
Jeotively  speaking,  a great  Improvement  him  been  obtained.  With  the 
uncompensated  system,  the  pilot  had  to  generate  an  acceleration  term 
in  the  crossover  region  in  order  to  satisfy  the  Crossover  Modal.  That 
is,  the  pilot  modo.1  had  to  become 

Vp  - Kp  e‘BT«  ■(■  + 1.111)  (1*1) 

which  is  very  difficult  for  the  human  pilot  and  loads  to  a high  pilot 
opinion  rating,  which  is  unfavorable,  and  a largo  pilot  time-dolay. 

The  compensator  allowed  the  pilot  to  adopt  the  form  of  a simple  first 
order  lead  with  a time  constant  of  from  one  to  two  seconds,  depending 
on  the  engine  time-delay.  This  rnsultn  in  a low  pilot  opinion  rating, 
which  is  an  indication  of  good  handling  qualities,  and  a small  pilot 
time-delay  (about  0.2  second). 

Pecommcndatlona 

An  experiment  using  the  analog  computer  and  a visual  display 
should  bo  set  up  so  that  the  control  design  can  be  tested  with  human 
pilots  in  the  loop.  This  would  check  the  accuracy  of  the  pilot  model 
used  and  would  reveal  what  improvement  is  possible  with  the  control 
systems  that  have  been  designed. 

Whan  attempting  to  model  the  aircraft  engines,  a need  was  found 
for  information  on  dynamic  modeling  of  aircraft  engines  of  various 
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types.  An  excellent  subject  for  a research  topic  would  be  the  dynamic 
modeling  of  both  turboprop  and  turbojet  engines. 

The  possibility  exists  that  the  ACLS  will  be  installed  on  air- 
craft with  more  than  two  engines  in  the  future.  A study  could  be 
performed  on  the  bunt,  method  of  differential  thrust  control  for  those 
type  of  aircraft.  For  example,  it  may  be  better  to  use  only  the  two 
outboard  engines  or  the  two  inboard  engines  rather  than  all  of  them. 
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Hoot.  Locui?  Plots  for  Control  Hyntcmn  with  Tima- Lag 

An  ortfiontlnl  part  of  the  design  vnu  the  use  of  root  locus  plots 
Of  ront.ro  1 nyitemn  Witli  timo-lng.  Theim  plots  enn  ho  easily  oht.nl nod 
from  the  ROOTL  program  in  the  library  of  AFITPROQUAMB  on  the  CDC 
CYDER/ 7^  digital  computer  at  AP1T.  The  principle  behind  plotting 
the  root  louun  for  control  systema  with  time-lag  was  found  in  Rof,  4. 

If  0(b)  ie  the  open  loop  transfer  function  for  a control  system, 
then  the  root  locust  is  obtained  from  the  characteristic  equation, 
which  la 


<](■)  * -1  ( A-l ) 

Thin  complex  mpml ion  may  he  opl it  in  two  parts  so  that  the  magnitude 
of  0(h)  is  equal,  to  one  and  the  angle  of  'l(n)  la  equal  to  tlflO  degrees 
(for  the  principle  plaint).  The  root  locus  is  a plot  of  tho  angle 
condition  and  the  gain  values  on  the  locus  are  found  from  the  magni- 
tude condition. 

If  tho  open  loop  transfer  function  has  a time-lag  equal  to  T, 


then  the  characteristic  equation  in 

G(s)  e”^fl  ■ -1  (A-2) 

Now,  the  angle  condition  becomes 

^a(s)-wT  - ± n 1 , n-1,  2,  3 (a-3) 


oi?WTir-5i 


In  order  to  draw  the  root  Ioqub,  a family  of  root  loot  are  drawn 
for  varioun  atiglen  of  (1(b)  and  a second  family  is  drawn  on  the  same 
o-plane  of  varioun  angles  of  -wT.  The  loci  for  angles  of  0(a)  are 
difficult  to  conn t runt , but  the  loot  for  -uiT  angles  are  simply 
horJ'.'.ontul  linen,  c ronn  I ng  the  Jm  axis  at,  right  nngluii.  Pnlnta  are 
now  plotted  of  tho  intersection  of  the  two  families  of  loci  where 
the  interoocting  ,'oci  values  total  1100  degrees , When  the  points 
are  all  Joined,  the  result  is  the  root  looua  for  the  time-lag  system. 
Tho  gain  valuau  can  be  determined  from  the  magnitude  condition  aa 
in  the  oanti  without  time-lag,  The  magnitude  equation  beoomoa 

| a(a ) e”^8  J ■ 1 (A-4 ) 

When  a control  nyutem  has  a time-lag,  there  are  an  infinite 
number  of  roots  of  tho  characteristic  equation.  The  root  looue 
contains  an  Infinite  number  of  branchen  when  tho  o-plano  is  extended 
to  infinity.  This  reaulto  from  the  repetition  of  the  angle  loci 
of  -wT  in  tho  above  construction  procedure.  The  infinite  number  of 
branches  lie  in  horizontal  bandn  as  seen  on  the  root  loci  in  Appendix 
B,  When  the  time-lag  increases,  the  bands  become  closer  together. 
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Appaiullx  P 

Hoot  U.'cv\p  Plot  w Unyrt  I'q r yftp  Pc it an 


Fig.  B.2  Root  Looub  Plot  for  Yaw  Rat*  Feedback,  t * 0.I4  Saoond 
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YAW  RATE  FEEDBACK,  TRU  = t .0 
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Fig.  B.T  Root  Locus  Plot  for  Lead  Compensator,  x = O.k  Second 


Boot  Locus  Plot  for  Lead  Compensator,  t =0.8  Second 
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Appendix  C 

MIMIC  Simulation  Program 
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♦♦♦MIMIC  S OU  RC  £>L  A NGU  A GE  V POGRAM*  *♦ 
ILIST 

UUCOMPtNS  AT  t.'I)  SYSTEM 
CONSTANTS 

CONtOAMP, START) 
CROSSOVER  PARAMETERS 

PAR(WC»TAUA,TAUE) 

HEADING  DISTURBANCE"  MODEL 


PULSE  GENERATOR  ( RESETS  TM.  PULSE  AMPLITUDE  TO 
A NEW  GAUSSIAN  NUMBER  { V..RY  .2  SECOND 


YES 


FAKE  s SRG ( ST  ART ) 

DEL  * SR1(T) 

Y'-S  ••  RSW  ( DEL,  FALSE*  TRUE,  FALSE) 

PULSE  * RUG ( 0. » 23. i 0) 


. FILTER  TO  S I MUL ATE  WIND  CU’ OFF 
« AT  1 RAD/ S:  C ON  POWER  fP.  UTRAL  DENSITY 


t 


WIND 

AIRCRAFT  TRANE 

01  ST  D 
OIST 


* INI  (1.54* ("ULSr-WIHD) , 3»> 
r L •'  FUNCTION  - HI* ADING/WIND 

* lNT(DAMp»v;nn-,07e.8*0ISTi}- 

* INT ( OISTD » . ,) 


0:«56*OICT,0.) 


CROSSOVER  MODEL 


ERROR  » r'IST-PSI 

ERROR  1 * NO*  I NT  ( r.  * » OK*  Q « ) 

TAU  « r ALIA  ♦ T A U t 

PSI  * 3R? (rRROF li TAU) 


AVERAGc  MEAN  SQUARcO  ERROR  ANALYSIS 


SERROP  * INT  (!  RROPV  RFOFfD.) 
MSERk  * GEPROR/T 


Bl 


GE/EE/T*t-51 


FINISH 


FIN(T, 10  0 • ) 

PRINT  INTERVAL 

OT  * .1 

INTEGRATION  STEP  SIZE 

OTMIN  * ,01 
OTMAX  « OTMIN 

OUTPUT 

MOR (TIME, PULSE, 0 1ST, PSI, ERROR, MSERR) 
OUT(T, PULSE, GIST, PS1 , ERROR , MSERR) 


PLOTTING  INFORMATION 


tPLOT  ONLINE ( SC l 

PLO(T,OIST,PSI) 

SC  A ( 1 * , . 0 J7  , , C 0 7 ) 
ZER(C.,5C . ,50.) 

TT  X ( T I ME  IN  SECONDS) 

TTY ( ANGLE  IN  RADIANS) 
OPT(1.40,1.  ,1.,) 


OE/ EE/74-51 


FUNCTION  SRI  ( A) 

* 

* FUNCTION  Sun PROGP AM  SETS  SRI  EOUAL  TO  ZERO 

♦ WHEN  SUBROUTINE  FIRST  CALLED  ON  A THE  EQUAL  TO 

* A MULTIPLE  OF  0.2  SECONDS 

♦ 

CO  MM ON/S WITCH /I  OUT, IPAR, DUMMY (1?) 

IF  (IPAR.  Nt:.l>  GO  TO  1 

Lb  0 

10  1 CONTINUE 

SR  1*  10. 

IF (A-C.  LT.Q. 0 01. AND. L.NE.O ) RETURN 
X*  A/  • 2+  , 0 00 1 
NX  “X 

15  Y*  X-NX 

IF  (Y.LE.  0.003  1)  $R1«0. 

1*1 
C*  A 
RETURN 

20  ENO 


FUNCTION  SR2(XIN, DELAY) 

* 

* FUNCTION  TO  DELAY  XIN  (TOTAL  T IMF  DELAY) 

* 

* 

DIMENSION  XSTORUS40) 

C O MM ON/ SWITCH/ IO UT , IPAR, OUMMY ( 17) 

INTEGER  OUT 

LOGICAL  TEST 

IF  (IPAR,  NE,  1)  GO  TO  101 

in*q 

OUTbO 

TEST", FALSE. 

DOIGO  1*1,640 

100  XSTOR(I)*0.0 

101  CONTINUE 

IF  UOUT.EO.  1)  TEST*. TRUE. 

IF  (Tf’.ST)  10,1 
1 NUM«4Q(1*  (DELAYS.  00001) 

SR2«0,Q 
RETURN 
10  CONTINUE 

IFdOUT.EO.l)  GO  TO  20 

IN*IN+1 

INI*  IN 

$R2*XSTOR  (INI ) 

0UT*NUM-IN1 
IF  (OUT. EQ . 0 ) I N* 0 
20  XSTORdNi  > *XIN 
RETURN 


UZ*/  ZiEi/  | 


FUNCTION  SP.5U) 


* FUNCTION  TU  INITIALIZE  GAUSSIAN  01  ST PlftOT ION 


DATA  N/0/ 

IF (U.Nl. 0 Hi  0 TO  1 
CALL  RAN3LTCA) 
iO  i N=i 

Sr>>0* 

RKTURN 

END 


B4 


J 


■ - — 
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♦♦♦MIMIC  SOURCE-LANGUAGE  Pi:OG PAM*** 

I LIST 

! COMPENSATED  SYSTEM 
! CONSTANTS 

CON ( CAMP*  ST  ART) 

. PILOT  PARAMETERS 

PAR  ( KSU9F, ? AUE , TSIIPL  , TSU3I ) 

. AIRCRAFT  AND  COMPENSATOR  PARAMETERS 

PAR(K,TAU\»Z1»P1»Z?»P2) 

. HEADING  DISTURBANCE  MOOEL 

. PULSE  GENERATOR  (RESETS  Tl'f‘  PUISL  AMPLITUDE  TO 
* A NEW  GAUSSIAN  NUM1ER  EV.RY  .2  SECOND 

FAKE  ■ 3RD ( START ) 

DEL  • SR  1 ( T ) 

YES  ■ FSW<OLL, FALSE, TRUE, FALSE) 

YES  , PULSE  • "’NO  ( 0 # , 2 J . t fl) 

*,  FILTER  TO  SIMULATE  WIND  ClHOFF 
, AT  1 RAO/ SEC  ON  POWER  S^-CTRAL  OENSITY 

HINO  ■ INT ( 1 • 54* ( pULSF» WIND ) i 0 • ) 

. AIRCRAFT  TRANSFER  FUNCTION  - HEAQING/WIND 

OISTP  * INT(r)AMP*WINP-.  C76A*MSTn-.  J?!5G*nTS 
OIST  ■ INT ( CISTD » T , ) 

! ACL S INNER  LOOP  MODEL 

C 

. PILOT  03SERVCS  HEADING  AND  DISTURBANCE 
PILOTI  ■ JIST-PSI 
*.  CROSSOVER  PILOT  MODEL  <YP) 

. K1UBP  IS  THE  PILOT  GAIN 

PI  LI  ■ KSU3P*PIL0TI 
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cie/ee/tMi 


. PILOT  DYNAMICS 

.(IF  TSUBI  IS  0..THEN  PlL2«'rUnL*DER (T ,PIL1, 0.  > +PIL1 
LAG  » 1 • /T  SUB  I 

HOLAG  * FSW(  TSUi3I  , FALSE  ,TRUE  , FALSE) 

YeSLAC.  * HOT  ( MOL  AG ) 

NOLAG  PIL2  ■ TSU3L*ltK'(T,PILl,0.) +PIL1 
YESLAG  P1L2  ■ L AG* ( TSUP L* P1L1 ♦ INT ( FIL t-PIL 2 , 0. )) 


PILOT  TIME  DELAY 

PILOTO  ■ 5R2 ( PXL2 1 TAUE) 

YAW  RATE  FEEDBACK 

PL Ai  ■ PILOTO-PSID 

COMPENSATOR 

COMPENSATION  CONTROL  LOOP  SCNSXTXVXT Y«< 

PLA2  * K*PLAl 

LEAD  COMPENSATOR  S+Z1/S+P1 

PLAT  ■ PLA?*INT < 71 *PL A 2-P 1* PL  A 3, n . ) 

LEAD  COMPENSATOR  S*Z2/S*P2 

PLA  ■ PLA3»INT<Z2*PLA3-P2*PLA|i).l 

LEAD  COMPENSATOR. < S f W 1 11 1 W (S* 1 1. Ill)  CHANGES 
ENGINE  DENOMINATOR  FROM  (5  + 1. 1111)  TO  (S+tl.lll> 

ENGINE  DYNAMICS 

NT  ■ TNT  ( PLA- i l-  lli*NT,0.  ) 

ENGINE  TIME  DELAY 

N « SR3 ( NT » T AUA ) 

AIRCRAFT  DYNAMICS 

PSIO  * INT ( N * 0 • ) 

PSI  ■ INT  (PSIOi  0.  ) 


AVERAGE  MEAN  SQUARED  ERROR  ANALYSIS 

ERROR  ■ INT ( PILOT  I* PZL OTZ»  0.1 
MSERR  ■ iiRROR/T 
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FINISH 

FIN(T*13C«) 

PRINT  INTERVAL 

OT  ■ .1 

INTEGRATION  STEP  SIZE 

DTMIN  • .01 
DTMAX  ■ OTHIN 

OUTPUT 

HOR{TIME»PULSF»DIST*PSIipSIO|HSEPP) 
OUT (T. PUL SE i DIS T uPSI fPLiD* MSCRRI 


PLOTTING  INFORMATION 


IPLOT  ONLINE (SC) 

PLO(T«niST*PSI) 

ZER(0. .SO.,30.) 

SC  A ( 1 * f ♦ C C?  » . 0 0'/  ) 

TTX ( TIME  IN  SEOONOS) 

TT  Y (ANGLE  IN  RADIANS) 
OPT ( 1 • A * 1 • »!•»?»  1 1 • ♦ 1 • ) 


. END 


END 
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FUNCTION  SRI (ft) 

FUHCTIOM  SUBPROGRAM  SITS  iPl  FOtJAL  TO  Z'R3 
WHEN  SU'ROUTInF  FIRTT  CALLAO  ON  \ TIMf.  EOJAL  TO 
A MULTIPLE  OF  0.2  SECONDS 

COMMON/ SWITCH /TOUT, T°ftR, DUMMY < \7) 

IT U PAP.Nt.l) 00  TO  1 

L*0 

10  1 CONTINUE 

SRl-iC. 

IF  (A-C.LT.O.  0 01.  ANO.L  .NO*  RETURN 

X*  A/ « 2+ . 0 00 1 

NX*X 

15  Y*X-NX 

XFIY.LE.  0,0  00  1)  SRi«0. 

L ■ 1 
C*  A 

RETURN 

20  END 


/ 

5 


FUNCTION  SR2(XtN» DELAY) 


5 


FUNCTION  TO  DELAY  XIN  (PILOT  TIME  DELAY) 


10 


19 


20 


29 


C 


30 


OlMf.HSION  XSTQR  (1J5) 

C CIMMON/CW  ITCH  / IOUT  * IPAR,  DUMMY  (17) 

INTt’GCR  OUT 

LOOT  CAL  TEST 

IF  (I  PAR*  NL.  1)  GO  TO  101 

ZN>0 

OUT*  0 

TEST*. FALSE. 

OOl':  0 I«  1 tl  35 

100  X $ T ( i R ( I ) ■ 0 . Q 

101  CONTINUE 

IF  (1  OUT.  EO,  1)  TEST  * .TRUE. 

IF  (TEST) 10,1 
1 NUH-400* (DELAY*.  O'  001) 

SR 2*  C. 0 
RETURN 
10  CONTINUE 

IF  (IOUT . EQ. 1)  GO  TO  20 

IN*IN*1 

INI* IN 

SR28XST0R  (INI ) 

OUT*NUM- INI 
IF (OUT. FQ, C ) I N*0 
20  X5TOR  (XfU  >-XZN 
RETURN 
END 
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UE/EE/7*»-31 


1 


10 


IB 


20 


) 


FUNCTION  SRJfXIN, DELAY) 


* FUNCTION  TO  3 EL  AY  XIN  ICN'iINF:  TIM*  DELAY ) 


niwrNSION  XSTORmO) 

CO MMON/SW ITCH /I  OUT , I PAR, DUMMY (17) 

I NT'.'  Gl'R  OUT 

LOGICAL  TEST 

IF (I  PAR*  NE, 1) GO  TO  101 

XN*0 

OUT-O 

TEST*. FALSE. 

DO  10  0 1*1, A00 

100  XSTOR  (I ) ■ 0, 0 

101  CONTINUE 

IF  (I  OUT • EO«  1)  TEST*, TRUE. 

IF  ( TEST ) 10  1 1 
1 NUH-ACO* (DELAYS. 00001) 

SR3-C.0 
RETURN 
10  CONTINUF 

IF  (I  OUT *EO,l)  GO  TO  20 

IN*IN*1 

INI* IN 

SR  T*X5TOR  (INI ) 

OUT*  NUM- I N1 
IF  (OUT.  ro , C ) I N*  J 
20  XSTORIXNi  )»XIN 
RF  TURN 
END 


FUNCTION  SRB(A) 


FUNCTION  TO  INITIALIZ'.' 


i 


GAUSSIAN  IT  STR  IlHiTION 


DATA  N/0/ 

IF  (N  , tic,  0 ) GO  ro  1 
CALL  RANSET(A) 

J 1 N*i 

SR->»  0, 

RL  TURN 


Appendix  D 


Pul  bo  Representation  of  Hand  Limited  White  Noise 


Hand  limited  white  noise  hns  a power  spectral  density  which  is 
flat  up  to  the  cutoff  frequency  and  then  In  zero  at  higher  frequencies 
«u  in  Klg.  D.l  . When  the  mean  of  the  Oauimlan  ill  n I.  r I but  Ion  for  the 


Klg,  D.l.  Power  llpoctrs'l  Density  of  Hand  Mint  tad  White  Noise 


noise  1 u zero,  the  variance  run  \w  found  from  th#  autocorrelation 
R(0),  which  in  nqual  to  the  area  under  the  power  Hpectral  density  curve 
from  w " - *’  to  w “ + •*,  divided  by  ’.,*1) . The  following  equation  nxprsssas 


this  (Rnf . 17 i 3dfl ) i 

W(0)  * na 


B(iu  )dw 


(D-U 


where  a .In  the  variance  of  t.hw  Onus si an  diutr I button.  Th»  value  of 
B()  in  Pig.  D.l  can  Im  calculated  from  the  relationship  In  equation 
(U-i)  and  Is 

8 .Jfli  <D-»> 

a 
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,’inrpHilitffe' 


The  band  limited  white  noise  la  represented  by  a series  of  pulses, 
whose  amplitude  leu  tittuosian  distribution  with  mean  equal  to  sero, 
and  standard  deviation  equal  to  a , and  whoso  width  is  T.  This  pulse 
train  Is  shown  in  Fig.  D.P.  The  autocorrelation  of  then®  pulses 
appears  an  in  Fig.  0.3. 




— 

— 

-• 

" 

e-  T 

— — 

— 

Time 

Fig.  I).?  Paine  Train 


Fig,  0,3  Autocorrelation  of  a Pulse  Train 
The  power  spectral  density  of  a process  Is  equal  to  the  Fourier 
Transform  of  its  autocorrelation , This  Is  expressed  as 


H(w) 


e“Juf  R(x)dt 


(0-3) 
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1SIBI 


FfYnT.  ' 


OE/EB/7^-51 


Therefore,  the  power  spectral  density  of  the  pulse  train  is 


B{aO 


r t 
J -T 


“JuiT 


O^dT 


2tJpl’ 


2 in  !?»»>  wT 


u)T 


(D-k) 


which,  when  plotted,  looks  like  Pig.  D.><.  Tf  the  pulse  width  Is 
kept  small,  the  center  part  of  the  curve  will  be  kept  flat.  Hie 
pulse  train  can  the'*  be  passed  through  a first  order  filter  to  cut 
off  the  power  spectral  Annuity  uurve  beyond  the  flat  portion.  The 
roault  is  a clone  approximation  to  the  ideal  band  limited  white  noise. 
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